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confinement in this theory. We also show how one can incorporate the physics of the anomaly 
in this theory. Thus, our approach is a good pedagogical tool to explain how pions arise in real 
QCD. Another advantage of our approach is that we can easily design efficient cluster algorithms 
to compute a variety of quantities close to the chiral limit, thus allowing us to understand the low 
energy physics in a QCD-like setting from first principles. 
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1. Introduction 

One of the outstanding problems in lattice QCD is to compute low energy hadronic observ- 
ables, which are dominated by the physics of light quarks, with controlled errors. Although 
progress is being made in terms of fermion algorithms for lattice QCD, it will still be difficult 
to approach realistic quark masses in the near future. As a result, calculations are performed at un- 
physically large quark masses and then extrapolations to realistic quark masses are performed using 
chiral perturbation theory. However, for such an approach to be reliable, we must know the range 
over which the chiral expansion is valid. Unfortunately, this is not very well understood. In fact 
a systematic study of chiral perturbation theory as an effective field theory that describes a more 
fundamental lattice field theory has not been attempted in many situations. Thus, it is useful to 
study a simpler lattice field theory with the same symmetries of lattice QCD and make connections 
with chiral perturbation theory. This is the main motivation for our work. 

We study the physics of strongly coupled lattice QED with two flavors of staggered fermions 
which we refer to as our model throughout this manuscript. Although this model does not describe 
full QCD, it is an interesting model which has the same symmetries of two flavor QCD and was 
recently used to study the chiral phase transition [[Ij]. Note that even a t/(l) gauge theory exhibits 
confinement in the strong coupling limit, and the taste symmetry is irrelevant since it is maximally 
broken at strong couplings. Our model allows us to study the pion physics from a fundamental 
lattice field theory, very similar to QCD, so that we can understand the usefulness of chiral pertur- 
bation theory as an effective description of low energy physics. It is important to note that unless we 
find a way to fine tune our model we will be dominated by lattice artifacts since at strong coupling 
the pion decay constant F n is naturally close to the cutoff. In order to circumvent this problem, 
we define our model in d + 1 dimensions where d = 4 is the space time dimensions. The extra 
dimension plays the role of a fictitious temperature which allows us to tune to a critical point where 
Fjt will be much smaller than the cutoff. Thus, we can still explore the physics of a continuum limit 
even in the strong coupling limit. 

The motivation for studying a strongly coupled theory is that the gauge dependent degrees of 
freedom can be integrated over, which significantly simplifies the theory. Further, a new class of 
algorithms, the Directed Path Algorithm, has recently been designed for strongly coupled gauge 
theories, by which it is possible to study the chiral limit very efficiently We have extended this 
algorithm to our model and are currently testing it. Here we report on current progress. 

We also wish to compute via Liischer's method [Q] quantities that are measurable in scattering 
experiments. Specifically, we would like to measure scattering lengths, characterize and understand 
resonances which may exist, and provide a setting to understand non-perturbative features that can 
arise in a QCD like theory and that go beyond the scope of chiral perturbation theory. Such a 
study appears impossible with the current technology of lattice QCD since some physical processes 
which occur for sufficiently light quarks, e.g. the decay of the p meson into two pions, is forbidden 
at the quark masses currently used. 

2. Model and Symmetries 

The Euclidean space action of the Nf = 2 QED model we consider is given by (Note that the 
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usual factors A are absorbed into the field definitions): 



'¥ x ¥ 



x+fl 



¥x+(l¥x 



mWxVx+^l ¥ x ¥x 



(2.1) 



where x denotes a lattice site on a d + 1 dimensional hypercubic lattice L, x L d . lff x and \}f x are 
two component Grassman fields that represent the two quark flavors of mass m, and <^, A is the 
U(l) gauge field through which the fields interact. Note that ji. runs over the temporal and spatial 
directions 0, 1,2, ...,d with denoting the temporal direction. The usual staggered fermion phase 
factors ri^x obey the relations: t]q x = T and T]f = 1 for i = 1,2,..., d. The parameter T is the 
fictitious temperature which will be used to the continuum limit. The coupling c will set the strength 
of the anomaly. 

We now discuss how our model has the same symmetries and symmetry breaking patterns of 
full QCD. It is first useful to note that any sum over lattice sites can be decomposed into a sum 
over even and odd sites. At c,m = 0, the action exhibits a global SU (2) x SU (2) x Ua(1) sym- 
metry. In particular, the action is invariant under the following Ua(1) and SUi{2) transformations 
(respectively): 
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SUr{2) is obtained by Vz, Vr and o 4^ Here Vj, and Vr are S?7 (2) matrices and can be param- 
eterized by: exp(/0 • a) where a,- is a Pauli matrix that acts on the flavor space. At c 7^ 0, Ua(1) 
is explicitly broken and the action is invariant under SUl(2) x SUr(2) x Z 2 . Thus, the coupling c 
induces the effects of the anomaly. Further at m ^ 0, it is necessary to set Vl = Vr for the action 
to remain invariant. Thus, with a mass term the chiral symmetry SUl(2) x SUr{2) is explicitly 
broken down to SUy(2). In order to mimic QCD we need to set c ^ and Hence, our model 

has the same chiral symmetry as full QCD. Further, based on previous mean field strong coupling 
calculations [BJ], we expect that the symmetry breaking pattern is also similar to full QCD. 



3. Mapping to a Monomer - Dimer - Pion Loop - Instanton Model 

The partition function of our model is equivalent to that of a classical statistical mechanics 
model involving configurations made up of gauge invariant objects such as monomers, dimers, pion 
loops and instantons We denote these as MDPI configurations. Note that a double monomer 
on a site breaks the C/a(1) symmetry but not the SUl{2) x SUr(2) symmetry and hence is called 
an instanton. In addition to these, each configuration can contain open loops of dimers which 
terminate on monomers, closed loop of dimers and oriented closed pion loops made up of oriented 
dimers. Explicitly, the partition function is given by: 

Z= £ Y\m nAx) m n " (x) c I{x) (3.1) 

[I,n' i ,n<',7$,n$,7t] l ] x ^ 
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where [I ,n ,n u ,71^,71^,71^] denotes a MDPI configuration. Note I(x) is the number of instantons 

% d 



on a site x, nj(x) the number of d monomers, n u {x) the number of u monomers, jrfj the number of 



d dimers, 71^ the number of u dimers, and 71^ the number of oriented ud or du dimers. The allowed 
values are: 

I{x)=0,2, n d (x)=0,l n u (x)=Q,l n£(x)=0,l fl*(x)=0,l Tufa) = -1,0,1 (3.2) 
Due to the Grassmann nature of the observables the following constraints must also be satisfied: 

J>J(*) = (3.3a) 



A 1 



*«w+E 



•n d (x) = 



Figure [I] gives an illustration of an MDPI configuration in 1 + 1 dimensions. 
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Figure 1: An example of a 2 x 2 lattice configuration as discussed in the text. 



4. Algorithm and Observables 

A Directed Path Algorithm can be constructed to update the MDPI configurations. Our algo- 
rithm is an extension of the algorithms discussed in [0, 01. We have three update routines: 



4 



Pion physics in two flavor strong coupling lattice QED 



D. J. Cecile 



1. A loop flip is an operation in which the orientation of a pion loop or a dimer loop is changed. 
Note that each dimer loop also contains two states. 

2. A loop swap which flips a dimer loop to a pion loop and vice versa. 

3. An update which creates and destroys instantons, double bonds, dimer loops, and pion loops. 
This is similar to the directed loop update constructed in [§, ||] and will be discussed in detail 
elsewhere. 

Numerous observables can be measured with this algorithm. The simplest are the three helicity 
moduli or current susceptibilities. In particular, for a conserved current /' (x), the helicity modulus 
(current susceptibility) is defined as: 




(4.1) 



where we are assuming a L t x L lattice. There are three conserved currents in our model. They 
are the axial, chiral, and vector currents which are given by: 

+ + (4.2a) 
(-lrWto-aJto] (4.2b) 
nfa) (4.2c) 

We also can measure correlation functions of pions defined as: 

(_iy+y _ 

G x (x,y) = — - — (W x e 3 Y.xW y o 3 Yy) ( 4 -3a) 
(-iy+y _ _ 

Gr,(x,y) = ^ — VMAV'vVA (4.3b) 
The corresponding susceptibilities, % n and Xt] are given by: 

* x,y 

The directed path algorithm allows a straightforward measurement of G(x,y) and %. Details will 
be given elsewhere. 

5. Results 

We have compared the measured observables to exact values on 2 x 2 lattices and found ex- 
cellent agreement as shown in Tables 1,2. Below are the analytic expressions for the partition 
function, the chiral, vector, and axial helicity moduli and the two susceptibilities, defined in the 
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T 


c 


m 


Algo. 


Exact 


Algo. 


Exact 


Algo. 


Exact 








yA 

w 




Y c 

W 




Y v 

W 




1.0 


0.5 


0.0 


0.8023(9) 


0.80246... 


0.5763(6) 


0.57721... 


0.5771(8) 


0.57721... 


1.5 


0.5 


0.0 


0.5212(7) 


0.52141... 


0.3275(6) 


0.32790... 


0.3274(7) 


0.32790... 


1.0 


1.0 


0.0 


0.7449(8) 


0.74534... 


0.5470(7) 


0.54658... 


0.5468(8) 


0.54658.. 


1.5 


1.0 


0.0 


0.4967(7) 


0.49720... 


0.3178(5) 


0.31821... 


0.3175(6) 


0.31821... 


1.0 


1.5 


0.0 


0.6667(8) 


0.66645... 


0.5016(5) 


0.50240... 


0.5014(6) 


0.50240... 


1.5 


1.5 


0.0 


0.4607(6) 


0.46147... 


0.3029(5) 


0.30325... 


0.3033(6) 


0.30325... 


1.0 


2.0 


0.0 


0.5812(7) 


0.58064... 


0.4519(6) 


0.45161... 


0.4521(6) 


0.45161... 


1.5 


2.0 


0.0 


0.4199(5) 


0.41927... 


0.2853(5) 


0.28451... 


0.2853(5) 


0.28451... 



Table 1: Helicity moduli for a 2 x 2 lattice as discussed in the text. 
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m 


Algo. 


Exact 


Algo. 


Exact 


1.0 


0.5 


0.0 


0.3601(3) 


0.359877... 


0.2172(2) 


0.217334... 


1.5 


0.5 


0.0 


0.2676(2) 


0.267764... 


0.1824(1) 


0.182429... 


1.0 


1.0 


0.0 


0.4040(3) 


0.403727... 


0.1429(2) 


0.142857... 


1.5 


1.0 


0.0 


0.2981(2) 


0.298322... 


0.1367(1) 


0.136731... 


1.0 


1.5 


0.0 


0.4220(3) 


0.422301... 


0.0801(1) 


0.080103... 


1.5 


1.5 


0.0 


0.3172(3) 


0.317264... 


0.0948(1) 


0.094767... 


1.0 


2.0 


0.0 


0.4194(2) 


0.419355... 


0.0323(1) 


0.032258... 


1.5 


2.0 


0.0 


0.3251(3) 


0.324754... 


0.0588(1) 


0.058961... 



Table 2: Susceptibilities for a 2 x 2 lattice as discussed in the text. 



previous section, for a 2 x 2 lattice (Note that c = c + m 2 ): 

Z(T,c,m) = 36T 4 + 64T 2 + 36 + c 4 + 12(1 + T 2 )c 2 + 8(1 + T)c 2 m 2 (5.1a) 
+ 327W + 16(1 + T 2 )m 4 + 48(1 + T 3 )m 2 + 32(r + T 2 )m 2 

2ZxY$(T,c,m) = 96 + 64r 2 + 16c 2 + 112m 2 + 32m 4 + 8c 2 m 2 + 32Tcm 2 (5.1b) 
+ 32(2r + T 2 )m 2 

2ZxY^(T,c,m) = 96 + 64r 2 + 16c 2 + 64m 2 (5.1c) 

64j2 2 

Y*(T,c,m) = Y$(T,c,m) + (5. Id) 

2Z x %Tt = 24r 3 + 16r 2 + \6T + 24 + 4( 1 + T)c 2 + 4(3 + AT + 3T 2 )c + 2c 3 (5. le) 

+ 8(1 + T)cm 2 + 16(1 + T + T 2 )m 2 
2Z x Xt] = 24r 3 + 16r 2 + \6T + 24 + 4(1 + T + m 2 )c 2 - 4(3- AT + 3T 2 )c (5.1f) 

- 2c 3 - 8(1 + T)cm 2 + 8(5 + 6T + 5r 2 )m 2 + 16(1 + T)m A 

The exact analytic results are shown in Tables 1,2 along with values calculated from our algorithm 
for various T,c and m. 
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6. Conclusions and Future Work 

Based on the results in the previous section, we conclude that our algorithm can in principle 
be used to study Nf = 2 lattice QED. Based on previous experience we expect the algorithm to be 
efficient for a range of parameters. We plan to study the model in d = 4 dimensions and match our 
results to chiral perturbation theory. We will assume that the lattice cutoff is of the order of IGeV, 
and tune T and c so that our F n and are close to their physical values. We then plan to make 
connections of our data to chiral perturbation theory as we change m. Ultimately, we also plan to 
compute the effects of the quark mass on pion scattering by measuring the appropriate two and four 
point correlation functions and extracting scattering phase shifts and lengths via Luscher's method 
as done in [|7|]. 
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